Anyonic Loops in Three Dimensional Spin liquid and Chiral Spin Liquid 
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We established a large class of exactly soluble spin liquids and chiral spin liquids on three di- 
mensional helix lattices by introducing Kitaev-type's spin coupling. In the chiral spin hquids, exact 
stable ground states with spontaneous breaking of the time reversal symmetry are found. The frac- 
tionalized loop excitations in both the spin and chiral spin liquids obey non-abelian statistics. We 
characterize this kind of statistics by non-abelian Berry phase and quantum algebra relation. The 
topological correlation of loops is independent of local order parameter and it measures the intrinsic 
global quantum entanglement of degenerate ground states. 
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I. INTRODUCTION 

Topological order of condensed matter system is the 
cornerstone of topological quantum computationfT]. The 
topologically ordered system have degenerated ground 
states which depend on the topology instead of symme- 
try. This degeneracy of ground state is robust against 
any local perturbation due to the gap between ground 
state and excited states 0- An interesting application of 
topological stablizer quantum codes is implemented in 
Kitaev's toric code modelQ. A topological color codes 
with quantum error-correcting capabilities is developed 
later [i I, this model removes the need for selective ad- 
dressing. 

There are different ways to perform quantum compu- 
tation in use of topologically ordered states. Usually the 
quantum gate or unitary transformation are implemented 
by brading quasiparticles. The transversal implementa- 
tion of quantum gate provides us another wayof topolog- 
ical quantum computation without braiding Q, in which 
the stablizer color codes is achieved in a 3-dimensional 
lattice model. While we focus on the braiding of loop 
quasiparticles in 3-dimcnsional spin liquid and chiral spin 
liuid in this paper. 

We established the 3-dimenisonal spin liquids by gen- 
eralizing Kitaev's two dimensional spin model. Kitaev's 
spin model on honjycomb lattice has a highly degen- 
erate ground state [gj. Abelian anyons and non-abelian 
anyons emerge in the background of short range spin 
liquid. These anyonic excitations obey exotic statistics. 
Quantum computation is based on anyon braiding. Dur- 
ing the braiding operation, anyons interact with one an- 
other only by encircling a topological nontrivial cycle. A 
proposal of the experimental operation of anyon braid- 
ing in Kitaev honeycomb model has been provided by use 
of ultracold atoms trapped in optical lattices 0], however 
there is some controversal on the details of this experi- 
mental proposal^]. 

Another most promising candidate of non-abelian 



'Electronic address: ' sity@itp.ac.cn| 



anyons may be realized in fractional quantum Hall states 
with V = 5/2[^[l3|. The fractional quantum Hall states 
break parity and time-reversal symmetry since fractional 
quantum Hall occurs in the presence of a strong magnetic 
field. Quantum information is stored in different fusion 
channels of non-abelian anyons [Tl. For two a quasiparti- 
cles following fusion rule a ■ a ^ 1 + ip, if fusion outcome 
is 1, we define the qubit as state |0). When fusion out- 
come is the state of the qubit is |1). We may use 2n 
antidots at which quasiholes are pinned to construct n 
qubits. 

Non-abelian anyons also arise from other two dimen- 
sional systems, such as Sr2RuO/^. In this p + ip super- 
conductor, the half-quantum vortices with flux hc/2e ex- 
hibit non-abelian statistics [ll|. Spin-1/2 chiral spin liq- 
uid states proposed by Kalmeyer and Laughlin possess 
fractionalized excitations[T5|. This kind of non-abelian 
anyons appears in 2D solvable spin models which are gen- 
eralization of Kitaev honeycomb model[l3|[3- These 
anyons are point quasiparticles confined in two dimen- 
sions. The main disadvantage of two dimensional system 
is that the expectation value of topological symmetry op- 
erators vanishes at finite temperature. This makes the 
topological memory thermally fragile 1^ 1^ . 

Since most physical materials in nature are three di- 
mensional, anyonic excitations most likely exist in three 
dimensions. One clue comes from a three dimensional 
exactly solvable spin model constructed by generalizing 
Kitaev's toric code to cubic lattice[l3|, in which the topo- 
logical order of its ground state is characterized by string 
net condensation as well as membrane condensation. The 
string net and membrane condensation also appear in 
another exactly solvable spin model [3, which is con- 
structed on a general 3-complex embedded in a closed 
and connected base manifold in three dimension. Its 
ground state degeneracy depends on the homology of the 
3-manifold. 

To establish three dimensional chiral spin liquid, we 
start from a different way to search for anyons in three 
dimensions by generalizing Kitaev honeycomb lattice 
model. Our interest focuses mainly on nontrivial loop 
excitations. Loops confined in two dimensions have no 
non-trivial entanglement without cutting each other. It 
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The paper is organized as follows: in section II, we con- 
structed exactly solvable three dimensional spin models 
on helix-decorated lattice. Their ground states are in- 
variant under time reversal transformation. A topologi- 
cal phase transition from gapless phase to gapped phase 
is observed. In section III, we established three dimen- 
sional chiral spin liquids by doping spins in the helix lat- 
tice. Their exact ground states break the time rever- 
sal symmetry. In section IV, the topological correlation 
among loops is studied from the point of view of topo- 
logical quantum field theory and Jones polynomial. In 
section V, non-abelian Berry phase of loops is presented. 
We studied the quantum algebra of loop statistics. 



FIG. 1: (a) The honeycomb helix lattice, (b) A helix spin 
chain, (c) The projection to X — Y plane. There are two 
types of elementary plaquettes, plaquette A and plaquette B 
. Plaquette A is the shadowed hexagon. Plaquette B covers an 
area of double hexagon, the oriented circles in the blank dou- 
ble hexagon denote the helix chains with different chirality. 
(d) The bonding configuration of the hexagon plaquette. 



II. EXACT SPIN LIQUID ON THREE 
DIMENSIONAL HELIX LATTICE 

We study the exactly soluble spin models. The exact 
ground state is short-range spin liquid. They all contain 
helix structure, so we classify them by the typical helix 
spin chain embedded in the three dimensional lattice. 



is in three dimensions that loops demonstrate much in- 
teresting statistical phenomena without breaking topo- 
logical constraint. There is a loop braiding group for 
the statistics of unknotted, unlinked closed stringsfioj. 
For links consisting of many entangled knots, their non- 
abelian statistics is closely related to the monodromy ma- 
trix in conformal field theory [2O, HH. Quantum loop 
gas applied to topological quantum computation is un- 
der rapid developing [22l|. 

There were some theoretical investigations on three 
dimensional chiral spin liquid based on effective field 
theory ^23^. However, no exact three dimensional chiral 
spin liquid is found so far to the authors' knowledge. We 
established a large class of exactly solvable three dimen- 
sional spin liquid and chiral spin liquid in this paper. In 
these models, there is only short-range interaction be- 
tween neighboring spins. The stable ground states of 
chiral spin liquid spontaneously break the time reversal 
symmetry. There exist loop excitations which may serve 
as basic qubits to store quantum information. Loop con- 
densation occurs naturally at ground state. 

The projection of these three dimensional spin mod- 
els to X — Y plane leads to soluble two dimensional 
models which are generalization of Kitaev honey- 
comb model. This projection operation does not keep 
time reversal symmetry. For example, the exact spin 
liquid we established on 'Cu'-sublattice of crystal green 
dioptase, which is the crystal structure of a material 
CuqSiqOs • 6i?20[23|, projects out the very two dimen- 
sional chiral spin model constructed in Ref. [l3|. 

Mandal and Surendran obtained a different exactly 
solvable Kitaev model in three dimensions (25j when we 
finished the research on the spin models in this paper. 
None of the exact spin models in this paper is identical 
to Mandal and Surendran's model. 



A. The Honeycomb Helix lattice model 

The first spin liquid model we proposed is established 
on coupled hexagonal helix lattice whose projection to 
the X — Y plane is Kitaev model on two dimensional 
Honeycomb lattice(Fig. 1). Two neighboring hexagonal 
helices have opposite chirality, three pairs of helix string 
surround one hexagon plaquette. The translational in- 
variant Hamiltonian reads 

i 

I T Z Z 

^ -'y"i-e6+2e2+e3"i-e6+e2-|-e3 
I J z z 

+ '^a;0'i_2e6 + 2e2+e3'''i-e6 + 2e2+e3}j (1) 

where the sum runs over i = mea + ne^ + /Gc with ba- 
sis vector Ba = (3V3a,0, 0), e?, — {^^a,^a, ^c) and 
Gc = (0,0, c)(Fig. l).(a and c are two different lattice 
spacings.) The Hamiltonian describes three dimensional 
coupled hexagonal helix chains. Each one chiral hexago- 
nal chain is equivalent to an exactly soluble one dimen- 
sional alternative coupling spin chain 

He = ^ JxO'2j'^2j+l + Jy'^2j + l'^2j+2 + Jz<y2j+2'^2j+3-i'^) 
2j 

The chirality of helix chain embeds in the operator prod- 
uct of the three component of spin. The chiral opera- 
tor a^a^a^ marks the right-hand helix and cr^cr^cr^ char- 
acterize the left hand helix(Fig. 1). Usually the non- 
vanishing chiral order parameter for spin liquid state is 
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defined as 

Oi]k = ^^i ■ i^i ^ CTk). (3) 

Another representation of the chir al op erator in terms of 
the conventional spin-1/2 fermions|26| is 

Oijk. — tijtjktki ^ik^kj^ji^ (4) 

where Uj = c|^Cj> as the hopping operator between two 
lattice sites. The chirality of helix chain may be charac- 
terized by the eigenvalues of 0^jk^ OijfejV') = 

The Hamiltonian commutes with two elementary loop 
variables. For loop A in Fig. 1(d), the conserved loop 
variable Pa is given by 

Pa - f^7f^8<^9f^?lf^i2'^?3- (5) 

And the conserved loop variable Pb in loop B (Fig 1(a)) 
is given by 

Pb = fTfa^fT^^CT^fTlfT^CTlCT^tTj^o, (6) 

which is equivalent to the product of two hexagon pla- 
quettes with their sharing bonds popping out as inter- 
layer bonding (Fig. 1 (a)). The eigen values of them are 
Pa = ±1 and Pb = ±1 because Pi = Pb = 1- No- 
tice that Pa and Pb may be simultaneously diagonalized 
because of [Pa,Pb] — 0. 

According to these conserved loop variables, we may 
obtain the exact ground state of this model by using Ma- 
jorana fermions to express spins [6], 

a? = z6?ci, fjf = zfefci, a? = ibfa (7) 

with constraint b^bybzC = 1. Majorana fermions 
bf,b^,bf, Ci obey the anti-commutations relation 

{bf,b^}^ 26^06,,, c2 = l, {bf,c,} = 0. (8) 

The nearest neighbor spin coupling term reads trf ct" = 
bfb'jCiCj, One may verify that the operator ufj = ibfb" 

commutes with the Hamiltonian and other u^- operators, 
i.e., 

K^.,H,]-0, [<^.,<] = 0. (9) 

Therefore, ufj may be viewed as the bond operator con- 
necting two nearest neighboring lattice sites and they 
may be though as 'good quantum numbers' which reflect 
the real physical conserved loop variables Pa and Pb- 
To see this point, we define the Wilson loop by a series 
of product of Uij along a closed path. There are two 
classes of elementary loops that can not be decomposed 
into smaller ones in this model, which are exactly the loop 
A and loop B in Fig. 1. The products of Uij along these 
loops are no more than the expressions of Pa and Pb 
in Majorana ferioms. Following Lieb's theorem(27|. the 
ground state must be vortex free state. The eigenvalue of 



the hexagon plaquette operator Pa at vortex free state is 
Pa = 1. This plaquette operator may be interpreted as 
the magnetic flux through the plaquette, so does the dou- 
ble hexagon plaquette operator Pb- The magnetic field 
line is perpendicular to the minimal surface of the Wilson 
loop. We may introduce magnetic helicity to character- 
ize the quantum entanglement between loops of different 
size. 



The Wilson loop operators are good quantum num- 
bers. They share the same Hilbert space with Hamilto- 
nian. Small loops do not gain extra energy when they 
expand into larger configuration. So these loops have no 
stress tension. They have the same ground state energy. 
Any spin configuration along a closed path which pre- 
serves Wc = — 1 for Wilson loop defines an elementary 
excitation. The excited state space is also highly degener- 
ate. There are two types of fundamental quasiparticles in 
honeycomb helix lattice model. Pa = — 1 and Pb = — 1. 
The honeycomb helix lattice model has nontrivial loop 
quasiparticles which obey exotic statistics in three di- 
mensions. Two loops confined in two dimensions could 
not entangle each other without cutting them up. But 
in three dimensions, they can avoid crossing each other 
when they wind around each other by keeping their topol- 
ogy. If they are linked with each other, the continuous 
transformation between different loop configurations cor- 
responds to the transition between different eigenstates. 
Quantum entanglement in the degenerate Hilbert space 
is the origin of non-abelian statistics. 



The the vortex free ground state is given by = 
-|-1. The Hamiltonian becomes quadratic Hh — 

hif^^juCi_fj,Cj,^^ where i,j are the site index of unit cell, 
/i, u indicates the internal degrees inside the unit cell. In 
the honeycomb helix model, the unit cell is a bounded 
pair of plaquette A and plaquette B{¥ig. 1). The funda- 
mental translational invariant vector of the unite cell is 
Vi = mea + neb + l^c with basis vector Ba — (3V3a, 0, 0), 
Gb = (^a, |a, ic) and Gc = (0,0,c)(Fig. 1). Perform- 
ing the Fourier transformations. 
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(10) 

we can derive the momentum representation of the 
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in the diagram covered by the inequahties: 
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FIG. 2: The phase diagram of the honeycomb helix lattice 
model. The shaded region is the gapless phase. 






FIG. 3: A cross section of the energy spectrum in — qy 
plane. The left panel is in gapless phase, the fermions fill 
in the low energy region while high energy region is cycled 
by low energy states. This means there exists closed Fermi 
surface. The middle panel is the critical region, the Fermi 
surface is nearly opened. The right panel is in gapped phase. 
The Fermi surface is completely opened. The dark area is the 
low energy region which is filled by fermions. The bright area 
indicates the higher energy region. 

Hamiltonian. The exact spectrum is 

EgiJ^, Jy, J^, q) = [Jl - 2 cos Z^/Zaq:cJlJl + 4 

+2cos(^y3aga; + - '^cq^)J^JyJz 

3 3 c 

+4cos(-\/3agj; + -aqy + -q^)J^JyJ^ 

3 /- 2,aqy cqz .,^47 
-Acos(-V 3aqx —)JxJyJz 

-2cos(^V3ag. - ^ + J,V. + SJj^ J| 

+4cos(cg,)jJj2j2 + 2cos(^V3a(7, - 

-2cos(^V3a,.-^)J,Vf 
+4cos(3agj, - {cq^) /2>)JxJyJt 

+2cos(3agy + ^)JxJyJt + J!)^ ■ (H) 



The ground state is highly degenerate due to innumer- 
able loop operators which are integral of motion. The 
whole Hilbert space are divided into two regimes — the 
gapped phase and gapless phase. The gapless phase falls 



I J.I' < 3| J, J, Jfcl + I J,f 



3| JxJyJz I 



\Jy? 



(12) 



where {i,j,k = x,y,z; i ^ j ^ k). The gapless phase 
region in the phase diagram is an infinite triangular area 
in the plane of |a;| — |j/|(Fig. 2). 

The quadratic Hamiltonian expressed by Major ana 
fermions operator has an equivalent representation in 
terms of complex fermion operators. We first express 
the spins by Majorana operators cr" — ib?Ci, and then 
introduce the complex fermion (30j on the cf coupling 
bonds, 



dr = (q + icj)/2, dl = {ci - iCj)/2. 



(13) 



The spin-spin coupling terms in Hamiltonian can be 
transformed into standard fermion representation. 



u%{2dld,~l), 

afaJ=4(4-Kdr)(4-e„-4-eJ 



(14) 



where r denotes the center of cfcrf bonds and u% — 
ibfb". The Hamiltonian ([1]) in the vortex free configura- 
tion is transformed into quadratic Hamiltonian expressed 
by standard fermions. The Fermi surface may be figured 
out from pT|) and has different topological behavior in 
gapped phase and gapless phase. We plot a cross sec- 
tion of the Fermi surface in Fig. 3. The center of the 
dark disc area is the lowest energy point. The bright 
area is the high energy region. The Fermi surface in the 
left panel of Fig. (3) is periodically arranged bubbles in 
three dimensional momentum space. Their cross section 
at qx — qy plane are closed circles without intersection. 
The Fermi surface in the gapless phase are topologically 
equivalent to compact sphere in momentum space. The 
middle panel is critical phase, at which the closed Fermi 
surface began touching each other but not connected. 
The right panel of Fig. 3 is in the gapped phase, the 
closed Fermi surface of gapless phase opened and con- 
nected each other (Fig. 3). 

The most natural topological index to distinguish the 
two phases is Hopf index which defines a map of from 
in momentum space to a sphere in coherent function 
space. The nontrivial two dimensional compact mani- 
fold embedded in the three dimensional momentum space 
comes from the gapless solution manifold E{p) — 0. The 
Hopf index is an integer which counts how many time 
the sphere in momentum space wraps the target sphere 
following the map fuopf ■ S^(j)) — *■ S'^. The nontrivial 
Hopf index means that the sphere enclose a monopole in 
momentum space. In the gapless phase, monopoles are 
enclosed by fermi sphere. They are at the right center of 
the dark disc areas in Fig. 3. They have either positive 
charge or negative. The monopole-anti-monopole pairs 
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FIG. 4: The second order derivative of the ground state spec- 
trum at Ja; = Jy = 1. In the region — 1 < Jz < 1 is gap less 
phase and outside is gapped phase. 



are fused into vacuum in the gapped phase, so we can 
not detect them through the Hopf index. The critical 
region occurs when the monopole-anti-monopole began 
touching each other but is not near enough to annihilate. 

The topological quantum phase transition from con- 
fined monopole phase to deconfined monopole phase 
leads to the divergence of the second order derivative of 
the ground state energy on the critical boundary. The di- 
vergence behavior of approaching the critical point from 
the gapped phase is totally different from that of gapless 
phase(Fig. 4). In the gapped phase, a sharp peak far 
from the critical point is observed in the region of Jz > 0. 
This peak is finite, it indicates a crossover transition. 

Now we turn to the spin correlations in the honeycomb 
helix lattice model. As shown in Ref. [2^ on Kitaev 
model of two dimensional honeycomb lattice, the two- 
spin correlation vanishes beyond nearest neighbor sepa- 
ration in different bonding directions. In the three di- 
mensional honeycomb helix lattice model, the two-spin 
correlation demonstrated the same behavior. The repre- 
sentation of spins by Majorana fermions factorized the 
eigenstates of physical Hilbert space as gauge sector and 
corresponding matter sector, |^) = |'0A/g)|V'g) with 

U,,\,P,)^brbf\^g)^u\^P,). (15) 

The spin operator af acting on any eigenstate adds a 
Majorana fermion at site i and one tt flux each to the 
plaquettes intersecting with the bond (ij)"- The two- 
spin correlation function 

w = {i^9\{i^M,wr{t)a^HM,m) (16) 

measures the probability we will detect the added Majo- 
rana fermion and flux at lattice site j. If Hilbert spaces 
of the two-spins at (i, j) have no overlap on the gauge sec- 
tor, the two-spin correlation is constantly zero. In this 
three dimensional honeycomb helix model, the two spins 
at (i,j) have no overlap in Hilbert space unless they are 
nearest neighbors. Therefore two-spin correlation van- 
ishes beyond nearest neighboring separation. 



The correlation between plaquette excitations is actu- 
ally multispin correlation. The non-vanishing multispin 
correlation comes from the operators which are product 
of an arbitrary number of the terms that appear in the 
Hamiltonian. We can compute the multispin correlation 
using the method in Ref. [28[ . Here we flrst consider the 
correlation between the two types of elementary plaque- 
tte excitations. As shown in Fig. 1, a spin starts from 
white lattice site i has three different propagating direc- 
tion. The interlayer spin correlation C"!^{t) are sharply 
cut off at nearest neighbor separation. This two-spin 
correlation measures the coupling between two Pb pla- 
quette excitations. The coupling between Pa plaquette 
excitation and Pb plaquette excitation is measured by 
three-spin correlation. One spin flip at i excites two pla- 
quette excitations. Pa = — 1 and Pb — —1. Their sharing 
boundary is along i j ^ k. For a spin flip propagat- 
ing along the path i ^ j ^ k, the nonzero three-spin 
correlation 

c^^'it) = (^.i(^M,|ar(t)af (i)a;:(i)iV'M,)i^,) (17) 

ends at another white site k. The two elementary pla- 
quette operators are the smallest loops, they can not be 
decomposed into smaller ones. As for the more com- 
plex loop-loop correlation, we shall introduce non-abelian 
Chern-Simons theory in the following sections. 

We introduce the other three dimensional helix models 
in the following. They have more or less similar property 
as the honeycomb helix model. We left the exact solution 
and present only the spin-spin coupling configuration. 

B. The Square Helix model 

The square helix model is composed of alternative XY 
spin chains each of which has four nearest neighboring 
chains coupled to one another through a^a^ bonding(Fig. 
5, (a)). The Hamiltonian includes two parts, Hsh = 
J2^xy + Hz, 

Hsh = Mj,k-i<^l,k + J-'^f3,k^fj,k+i] 

ij^k—even 

+ X! [■Jz^i~1.3,k'^l3,k] 

iin,j,k±2n 

+ X! [J^'^lj.k+l'^lj+l.k+l]- (18) 
i,j±n,k±2n 

A single helix chain is exactly the same single XY spin 
chain in two dimensional Kitaev model, which is equiv- 
alent to one dimensional Ising model with transverse 
field [1^. The multispin correlation function along a sin- 
gle helix chain demonstrates a long range order, 

= lim (H - [1 - iJy/J.n^'- (19) 

k=l 

In the three dimensional square-helix model, it is the 
product operator of the bonding terms in Hamiltonian 
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FIG. 5: (I)-(a) The square helix lattice model. (I)-(b) The 
triangle helix model on Cu-sublattice in Cu^SieO^ ■ 6H2O. 
(II)- (a) An equivalent demonstration of square helix model 
on cubic lattice. (Il)-(b) The deformed triangle helix model 
on cubic lattice. 



that has nonzero multispin correlation. While two-spin 
correlation in real lattice space vanishes at the nearest 
neighboring separation [2^. The chiral operator is the 
combination of different product of three spin operators, 
it is a good quantity to characterize the quantum entan- 
glement of three spins. 

The Hamiltonian (fTSl) is exactly soluble following a 
similar procedure as that in honeycomb helix model ex- 
cept a different gauge fixing. There are two types of 
fundamental plaquettes in square helix model. It is more 
clear to see the picture by deforming the square helix lat- 
tice into cubic lattice(Fig. 5). The energy minimum of 
ground state is achieved byO — tt — — tt-- - alternative 
flux pattern following Lieb's theorem. The integrals of 
motion Pa take a phase 27r for all irreducible plaquette 
A, Pa — erf CT|cr|crf cr|CT|crf crfcrf ct^'q — e^^'^ , and for all 

plaquette B, Pb = (^i(^2^t^t'^t^t'^7'^t — 

Elementary plaquette excitations are produced by flip- 
ping odd number of spins along the plaquettes A and B. 



At the excited states, the phase pattern of the two pal- 
quettes are Pa = and Pb = +1. The highly degen- 
erate eigenstates are classified by different spin configu- 
rations. The spin flip at lattice site i excites four tt flux 
each in its four neighboring plaquettes. Four plaquette 
excitations are always generated and annihilated at the 
same time. So there are many degree of freedom to oper- 
ate on them. Two of the four plaquette excitations may 
share a strong four-spin correlation in /c-direction. 



k+l.,k+2,k+'i 



I y z X X \ 

\"ij.k"ij,k+l"ij.k+2"ij,k+il 



(20) 



This four-spin correlation only covers four neighboring 
lattice sites. 



^fe,fc+l,fc+2,fc+m(^) — 0, m > 3. 



(21) 



While another two of the four plaquette excitation are 
connected only by nonzero two-spin correlation. Thus 
the correlation between plaquette excitations is different 
along different direction. 



C. The triangle helix model on Dioptast Lattice 

The triangle helix model is established on a three di- 
mensional lattice which exist in nature. The Cu atoms in 
dioptase CuqSIqO^ ■ 6H2O construct the exact triangle 
helix lattice, Cu-atoms form chiral chains along z-axes 
and project a honeycomb configuration to the X — Y 
plane [24 1 . 

The exact solution of the Heisenberg spin-spin cou- 
pling hamiltonian on the dioptase lattice is hardly reach- 
able. We derived an exactly soluble model by intro- 
ducing Kitaev-type spin-spin coupling interactions along 
different bonds on this three dimensional lattice(Fig. 
2(b)). The projection of the triangle helix model to 
X — Y plane is a two dimensional chiral spin liquid 
model constructed in Refs. [11 [11 which has a sta- 
ble ground state spontaneously breaking time reversal 
symmetry. While our three dimensional triangle he- 
lix model is invariant under time reversal transforma- 
tion. The two elementary plaquettes, octagon Pg and 
dodecagon P12 = erf cr|(T|cr4cr|fT|tT7tT|crg(TfoO'ficrf2, are 



both in TT— flux phase at ground states. Unlike the 
quantum phase transition from anti-ferro-magnetically 
ordered state to a quantum spin liquid in magnetic crys- 
tal CueSieOs ■ 6H2O, the quantum phase transition here 
in the triangle helix model is between two topologically 
distinguished spin liquid states. 



D. The Square-Hexagon-Helix Lattice Model 

The square-hexagon-helix lattice model projects one of 
the two dimensional models[l3| with Kitaev type's cou- 
pling in X — Y plane. The hexagon helix with opposite 
chirality are separated by isolated squares. The spatial 
distribution of the helixes is a pillar array on Kagome 
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lattice(Fig. 6). Further observation shows, the square- 
hexagon-hehx model is composed of three coupled Ising 
models which are established on multilayer Kagome lat- 
tice. 




FIG. 6: (a) The square- hexagon- helix Lattice model, (b) 
An equivalent configuration of square-hexagon-helix Lattice 
model on multi-layers Kagome lattice. 



III. THREE DIMENSIONAL CHIRAL SPIN 
LIQUID MODELS 

In the sections above, we have derived four typical spin 
models which are exactly soluble in three dimensions. 
Their ground states with short range interaction arc spin 
liquid states. It has been a long history to look for three 
dimensional chiral spin liquid, so far there is still no exact 
soluble three dimensional model to confirm its existence. 
In this section, we shall construct a series of three di- 
mensional chiral spin liquid models. They are exactly 
solvable. Their exact stable ground states spontaneously 
break the time reversal symmetry. These chiral spin liq- 
uids are likely to be realized by ultracold atoms trapped 
in cubic optical lattice. 



A. Chiral Square Helix Model 

We add two spins at all of the horizontal a^a'' bonds 
in the cubic version of deformed square helix lattice(Fig. 
5 (II)- (a)). One is at 1/3 of one lattice spacing and the 
other is at 2/3. We shift the added spins vertically a half 
lattice spacing along z-axis. If the spins are assumed to 
have only the nearest neighboring interaction, we get the 
first chiral spin liquid in three dimensions(Fig. 7 (a)). 
The Hamiltonian of the chiral square helix reads 



ilr.a — — 



ij ,k—even 



-even 

E[J„cr^ 1 . , cr^. , , 1 -f J^a^ 1 . , cr^. , i 

i±nj,k±2n 



,j,k'^i-^,j,k 



y 

I — 



i,j±n,/c±2n 

'^^'^1,3+ i ,fc+l + 1 ,fc+l 



J,, a 



ij+§,fc+l'^j.j + l,A;+f 



(22) 



Our interest focus mainly on the gauge symmetry of exact 
ground state. As shown in Fig. 7 (a), there are three el- 
ementary plaquettes: a triangle and two polygons. One 
of the two polygons has sixteen bonds, the other cov- 
ers twenty bonds. The energy is minimized by uniform 
flux configuration. The elementary plaquettes are all in 
TT-fiux phase at ground state. The tt— flux of triangu- 
lar plaquette flips under time reversal transformation. 
Although the two polygon plaquettes are invariant un- 
der time reversal transformation, the global flux pattern 
spontaneously breaks the time reversal symmetry. Thus 
the ground state is highly degenerate chiral spin liquid. 

The fractionalized elementary excitations are loop ex- 
citations in the three dimensional chiral spin liquid. 
Large loop operator is the product of the elementary pla- 
quettes covered by its minimal surface. If the loop has 
odd number of links, there must be odd number of tri- 
angular plaquettes confined in the area surrounded by 
the loop. A spin flip adds 7r-flux each to the plaque- 
ttes intersecting at that point. To some extend, the time 
reversal transformation behaves like transition operator 
between excited states and ground states. The loop with 
odd number of particles spontaneously breaks the time 
reversal symmetry. 

The correlation between plaquette excitations is mea- 
sured by multispin correlation which are spatially 
anisotropic. For example, a spin flip at site q excites four 
plaquettes vertices sharing a bond along z-axis. Two of 
them share a strong six-spin correlation 



<-'qqiq2q3q4P ~ \'^q^qi^q2^q3^q4^p/ 



(23) 



While another pair of plaquette excitations are only con- 
nected by two-spin correlation which are sharply cut off 
beyond nearest neighbors. This suggest that the fusion 
rules for plaquette excitations in this mode might be more 
complex than that in two dimensions. 
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(a) (b) 

FIG. 7: (a) The chiral square helix model on cubic lattice, (b) 
The chiral dodecagon helix model. The character structure is 
the dodecagonal helix. 




FIG. 8: (a) The chiral square helix model on the original 
square helix lattice, (b) The chiral spin liquid model on the 
triangular decorated dioptast lattice. 



B. The general principal to construct other 
cirtificial chiral spin liquid model 

There are many other ways to construct three dimen- 
sional chiral spin models with Kitacv's type of coupling, 
which requires that each lattice site branches out only 
three different bonds corresponding to OxC^x-, Oy^^y and 
OzOz- For example, we add two spins at every bond to 
divide it into three equal sections, and assume spins only 
interact with nearest neighbors, it naturally leads to a 
chiral spin liquid in three dimension. In Fig. 8, we 
present two examples, the triangular decorated square 



helix model and triangle helix model. If one deforms the 
triangular decorated square helix lattice model(Fig. 8 
(a)) into cubic lattice, the triangular decorated square 
helix is exactly the chiral square helix model. One can 
also construct other models by dividing some selected 
bonds into n(n = 2,3,4...) identical pieces and adding 
spins on the adjoint points with properly arranged bond- 
ing direction. The triangular decorated triangle-helix 
model(Fig. 8 (b)) is one of such constructions, to con- 
struct this model, we add one spin each on the two bonds 
branched from the same lattice site. 

Another different way of constructing is to spiral poly- 
gon plaquettes up. For example, we construct the do- 
decagon helix lattice model which has the same two di- 
mensional projection io X — Y plane as the diopaste lat- 
tice. Here it is the dodecagon that spirals up with trian- 
gles as bridge between dodecagonal helixes(Fig. 7). We 
can also obtain this chiral spin liquid model by replac- 
ing each lattice site in the honeycomb helix model(Fig. 
1) with a triangle. If we keep replacing the lattice sites 
with triangles time and time again, it finally leads to the 
fractal lattice model. This fractal model is different from 
Bethc lattice. There is no closed loop in Bethe lattice, 
but in this fractal lattice model, it is all sorts of loops 
that construct the whole lattice. These loops form the 
complete set of integrals of motion. Scale invariance lays 
in the heart of these exactly soluble fractal lattice models. 

Now we see that three dimensional chiral spin liquid 
emerges following a rather simple doping construction. 
The only thing we need is local interaction between near- 
est neighboring spins. When it comes to existing material 
in nature, it is possible get chiral spin liquid by doping 
spins in conventional three dimensional lattice. 

IV. TOPOLOGICAL CORRELATION OF 
ENTANGLED LOOPS IN THREE DIMENSIONAL 
CHIRAL SPIN LIQUID 

In the exactly soluble spin models obtained above, a 
closed path C consists of a series of bonds. Defining 
hjj = (yjiy'j for a given bond 7, the product of hj^ along 
C is a Z2 gauge invariant quantity which is nothing but 
the Wilson loop Wc = W^^c^lj- Hilbert space is 
divided into many gauge equivalent class {|V'c)} corre- 
sponding to Wc with 

Wc\^kc) = M^c). (24) 

Quantum information is encoded in this Hilbert space. 
The gauge invariant stats can only be distinguished by 
global operation, which protects the quantum informa- 
tion from decoherence due to local quantum entangle- 
ment with environment. 

Non-abelian excitations in three dimensional chiral 
spin liquid serve as elementary qubits to perform quan- 
tum computation. Basic logic gate may be constructed 
by designed braiding operation. To perform topological 
quantum computation, we shall operate on excitations 
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obeying non-abelian braiding statistics. Point quasipar- 
ticles in three dimensions are either fermion or boson and 
thus they are not good candidates. It is the loop exci- 
tations that satisfy non-commutable relation during the 
unitary evolution in ground state subspace. 

We first consider the unknotted loops which arc not 
entangled or linked with each other. If the loop corre- 
lation is extended to strong and weak coupling limit of 
lattice gauge theory, it shows that the loop correlations 
Sap '■= (Wc^Wcfi) decay exponentially with distance. 
At high temperature, the decay width is governed by the 
area A encapsulated by the two loops, Safi oc e'^^. At low 
temperature, the total perimeter \C'a \ + \Cp\ determines 
the decay behavior, S'q./j cx e'^^l'^^l+l'^'^l^. 

Loop-loop correlation is actually bipartite multi-spin 
correlation which can be calculated using the method 
developed in Ref. The area Aa enclosed by loop 

Cq. is composed of a collection of elementary plaquettes. 
The Wilson loop is just the product of all the elementary 
plaquette operators 

Wc^^YlPr, (25) 

i 

a indicates the type of elementary plaquette. For exam- 
ple, a = 1,2,3 in the chiral square helix model, Ci and 
C2 represent the two polygons and C3 represents the tri- 
angular plaquette. The degenerate Hilbert space of loop 
is the direct product of the subspace of all the elemen- 
tary plaquettes in Ac- For a given loop Ca, there exist 
many different Aa such that dAa = Cq. The dominate 
contribution to the correlation of two loops comes from 
the minimal surface M encapsulated between them. The 
boundary of the minimal surface is just the two loops, 
dM = Ca + Cp. We denote its Hilbert space as \1/m, the 
two loop correlation reads 

Sa^it) = {^i,\e^^"-^'>'WcAt)Wc,)\^a), (26) 
which has an explicit expression 

Sapit) = (vl/,|e^(^"^)* n n 1*^.)- (27) 

Expressing the spins in terms of Majorana fermions, we 
may divide the Hilbert space into different gauge sectors 
according to the static Z2 gauge field operators Uij = 
ib(bj. The multispin correlation can be computed on 
the gauge sectors [2^. However this computation is hard 
to tell the physical picture of topological entanglement 
between loops. 

In order to see the topological correlation more clearly, 
we take a more physical picture in the continuum limit. 
We view a loop as an electric current source radiat- 
ing magnetic field. The topological correlation between 
loops is equivalent to the topological entanglement be- 
tween magnetic field circles. Obviously no matter how 
far the loops are separated, the magnetic field lines are 
always entangled. When it comes to the knotted loops 



which are linked, entangled one anther, besides the con- 
ventional loop correlation which decays as their spatial 
separation increases, there exist topological correlation 
which relies on the global topology of links instead of 
the average distance between them. The topological cor- 
relation is immune to local fluctuation in space. Each 
knotted Wilson loop operator projects a subspace, this 
topological correlation is a good quantity to measure the 
quantum entanglement between different sectors of the 
total Hilbert space. 

To study the topological correlation of linked, knotted 
loops, we first take a preview in the continuum limit. 
Considering the Hopf link, in which two unknotted rings 
locked each other up, there is only one crossing point 
between loop Ca and the minimal surface of loop Cfj — 
dAp. The two loops have no intersection. For a more 
general case, one loop winds around another many times, 
the winding number is identical to intersecting points 

Car\ Afi = {P„i^ 1,2,. ..,m}. (28) 

The number of intersecting points is identified with a 
topological number called Linking number which can 
be calculated by Gauss link integral of lattice gauge po- 
tential. The three dimensional Jordan- Wigner transfor- 
mation mapping spin to fermions defined a gauge trans- 
formation, 

[/ = e*^J^'^°'^(j-'\ (29) 

The gauge potential is introduced as — V^tJ. The 
corresponding magnetic fieldjs^ is 

B{i,j,k)^V xAii,j,k) (30) 

with discrete derivative 

diA{i,j, k) = - A^^j^k- (31) 

The integral of the gauge field along Caa is Wilson loop 

The physicists' approach to the linking number of two 
loops Ca and Cp is using analogy with Ampere's law. 
Gauss link integral on the two entangled loops does not 
depend on the physical meaning of vector field tangent 
to them, they are valid for any physical vector field. So 
we think of loop Cq as being a wire carrying an electric 
current /. The line integral of the generated magnetic 
field B around the closed curve is AirLkl, where Lk is the 
linking number which count how many times the current 
winds around another closed loop. The linking number is 
just number of nontrivial intersecting points between Cq 
and the minimal surface of C/3 . One may also deduce this 
directly from the Biot-Savart Law. Therefore we define 
the topological correlation as 

S^pit) cx m gap cx Lk gap (32) 
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with Lfe as a topological index. A spin flip at the cross- 
ing point Pi excites the two loops at same time, thus 
they have overlap in the gauge sectors of Hilbert space. 
Then we derived a nontrivial contribution to the mul- 
tispin correlation function of two loops. This topolog- 
ical correlation survives local fluctuation and are valid 
for thermal states. It should be pointed out here that 
some intersecting points would disappear during topo- 
logical transformation. Those points do not contribute 
to the topological correlation. In the following when we 
mention intersecting points, it always means the non- 
trivial intersecting points which can not be eliminated 
by topological transformation. 

Knotted Wilson loops in the three dimensional spin 
models above also commute with Hamiltonian, they are 
good quantum numbers. The topological correlation be- 
tween linked knots can be described by topological quan- 
tum field theory. Witten found the quantum field theory 
to describe a wide class of topolo gica l invariant from the 
non-abelian Chern-Simons theory [20|. For a link contain- 
ing a collection of knotted, linked loops in three dimen- 
sional manifold L = {Ci,i = 1,2,3...}, we can compute 
the multiloop correlation function between its component 

Z{L) = (Q- • -C,) = Z{M)-^ I [VAy^ n WiiC), 

i=i 

(33) 

where Wc = rr(Pe'^ ^"''t") is the Wilson loop. S is 
the non-abelian Chern-Simons action S = ^ Jj^ Tr{A A 
cIA+^AaAaA), the first term in the non-abelian action 
represents the local magnetic Helicity 

Helicity = ^ijkBidjAk- (34) 

Z{L) is a gauge invariant which provides us a topological 
quantity to classify the geometric configuration of linked 
loops at different eigenstates. 

We compact the three dimensional lattice manifold 
into in the continuum limit. A sophisticate link em- 
bedded in this may be decomposed into two parts 
by cutting the using an as section manifold. We 
choose a proper cutting point so that the left piece Ml 
contains most complicated stuff of the link, and keep the 
right piece Mn as the simplest link configuration. The 
Feynmann path integral on Ml determines a vector \x) 
in the physical Hilbert space, the Feynmann path inte- 
gral on Mn determines a dual vector \ip). The partition 
function is the natural paring of the two vectors 

Z{L+) - (xlV), (35) 

L+ denote the over crossing of the two loops on the M^. 
The separated two parts of link puncture four points on 
the cross section manifold. The four punctured points 
represent two spin fields a and their corresponding an- 
tiparticles. Braiding these quasiparticles on leads two 
different configuration manifolds of the link, which we 



denote as Mi and M2. The corresponding path integral 
on Ml and M2 contributes another two different vectors 
\4>i) and |i/'2)- The two vectors can be viewed as be- 
ing generated by applying braid group operations on the 
original dual vector lip) of Mji. In topological quantum 
field theory, the Skein relation is encoded in 

a(xlV^)+/3(xlV'i)+7(xlV'2) = 0. (36) 

Namely, 

aZ{L+)+(3Z{LQ)+jZ{L-)=0. (37) 

Lq denotes the zero crossing in Mi and L_ denotes 
the under crossing in M2- Here the three parameters 
{q!,/3,7} depends on the special models. This relation is 
the fundamental constrain on the multi-loop correlation 
of sophisticate links in our exact three dimensional spin 
models. 

The topological correlation is invariant no matter how 
large one loop expands continuously. Any local change 
that break the topology of a link during evolution would 
sent out a signal by correlation function. The linking 
number Lk would jump from m to m ± i, (i 7^ 0), so does 
the topological correlation function 5^ cx Lk- If the loop 
breaks apart at the crossing lattice site, the topology 
of link changed, this change would be detected by the 
topological correlation function. 

As one loop expands without cutting the other loop, 
the eigenstates correspondingly evolute from a state of 
one loop to another in Hilbert space, but the topolog- 
ical correlation remains invariant. Therefore local per- 
turbation does not destroy the topological signal if the 
two loops does not intersect. This topological quantum 
entanglement does not depend on temperature, so it is a 
very good candidate for quantum memory. Three dimen- 
sional materials in natural are very promising for topo- 
logical quantum computation. 

V. STATISTICS OF STRINGS IN THREE 
DIMENSIONS 

The statistics of point particles is characterized by 
the phase change of many body wave function dur- 
ing the exchange of particles, for bosons, \4>{xi,X2)) = 
e'^^'^\ip{x2,xi)), and for fermions, \ip{xi,X2)) ~ 
e^'^\ip{x2,xi)). We are led back to the initial state 

\iPixi,X2)) ^e'^'<^\yj{x2,xi)). (38) 

After performing the exchange twice. The anyon statis- 
tics falls in the range of phase 7 varying continuously 
from 1/2 to 1. In three dimensions, any wrapping path 
of one particle around another can be topologically con- 
tracted to a point, that means two interchange is equiva- 
lent to an identical map, so the statistics of point particles 
is either fermionic or bosonic. 

The emergent string excitations bear exotic statistics 
in three dimensions. For the simplest illustration, if two 
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(a) (b) 



FIG. 9: (a) one loop wraps another loop, the chirality of the 
two loops in the upper is opposite and is same in the lower 
half- plane, (b) The exchange of two loops. The upper diagram 
depicts the world sheet of two loops with one loop flips its 
chirality during the exchange. The lower diagram depicts the 
case that two loops keep their chirality invariant during the 
exchange. 



unentangled and unknotted vortex rings leap through 
each other, a nontrivial statistical phase may emerge as 
Berry's phaseisT]. The leapfrogging vortex rings are a 
vivid demonstration of loop braid group [19j| for a collec- 
tion of miknottcd, unlinked loops. 

When it conies to discrete loops in three dimensional 
lattice, many exotic phenomena pop out. Any closed 
loop in the three dimensional exactly solvable models 
discussed above commute with the corresponding Hamil- 
tonian. The string operator along an oriented path Ca 
reads 

Oo - n (39) 

O denotes the counterclockwise order of operator prod- 
uct. Its complex conjugate string operator = 
n"gcO '^T has a clockwise product order. Since 

[Oa,H] = 0, the ground state of loop is an eigenstate 
of Oa following the equation 

OalV-a) =^l^a),2 = ±l. (40) 

In the degenerate Hilbert space of ground state, an 
initial state may end up with different final states after 
adiabatic evolution along a closed path. The conven- 
tional Berry phase needs to be generalized to a matrix 
of inner product between degenerate states. Notice that 
a loop sweeps out a two dimensional world sheet, the 
Berry phase (j) = {'ip{R)\\7 ii\ip{R)) of point particle along 
an one dimensional trajectory must be extended to cover 
the extra dimension of loop Cq.(s), where s measures 
the distance along the curve. We start from the gen- 
eral eigenfunction of time-dependent Schrodinger equa- 
tion, an adiabatic evolution in the degenerate subspace 



follows 

\MC{s))) = e-'/o"^(*)e»0=^|^^(c(s))), (41) 

where 

0a/3=i / dt [ dF(Va|9t^/3) (42) 
Jo JT 

is the extended Berry phase. We first consider an 
unknotted, unlinked loop Ca- The minimal surface 
Ea(i) bounded by Ca subtends a solid angle ila{R,t) 
at point R. A collection of solid angles of different loops 
{V,a{R,t),i}i3{R,t), ■ ■ ■ V,j{R,t)} constitute one geomet- 
ric representation of loops. In the continuum limit, the 
gradient of the solid angle field — V51(i?, t) induces an ef- 
fective magnetic field which circulates around the loop. 
This magnetic field could be computed from the Biot- 
Savart law. The extended Berry phase in Eq. has 
explicit expression in terms of solid angle field, 

0a/3 = J / dt [ dV{i:a\i'f3)dMR,t). (43) 

Jo JT 

If we assume the probability density pap = {4'a\'4'(3) is 
homogeneous, the quantum mechanical phase is propor- 
tional to the volume which is encapsulated by the world 
sheet bounded by the two loops Ca and C^. 

Exotic statistics phase is strongly dependent on the 
chirality of loops. For example, consider a red loop and 
a blue loop with the same orientation(Fig. 9), the ex- 
change of the two loops sweeps out a torus surface, the 
statistical phase is proportional to the volume of the torus 
^0/3*^"* • chirality of one loop flips during the ex- 

change, the world sheet becomes Klein bottle which is 
unoriented manifold with only one side. There is no def- 
inite volume integral for the quantum mechanical phase 
in this case. Here we take another equivalent but sim- 
pler exchange procedure by wrapping the solid red loop 
around the dashed blue loop. The dashed blue loop en- 
closed in the torus surface of world sheet behaves like 
electric current which produces a magnetic field follow- 
ing the Biot-Savart law 

dBb = CQr~^dsb X r. (44) 

The volume integral enclosed by the torus is computed 
from 

V = J^^ = JdBbXdSr- dR. (45) 

If the two loops have the same chirality, the statistical 
phase factor is proportional to the volume of a torus 
world sheet, it is positive. On the other hand, if the loops 
have opposite chirality, the volume integral is negative. 
Generally speaking, two loops with opposite chirality be- 
have like fermions, while two loops with same chirality 
obey bosonic statistics. 

Another key element governing the statistics of loops 
is the total number of particles confined along the loop. 
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The vortex excitation in the minimal surface surrounded 
by the loop behaves like quasiparticles of magnetic flux. 
We simplify the statistics of two loops as that between 
two quasi-particles. If the sum of spins of one loop is half 
integer, the quasiparticle is fcrmion. If the total spin is 
integer, the quasiparticle is boson. Moreover, quasipar- 
ticles must carry the orientation of loops, so they have 
internal degree of freedom which contributes nontrivial 
statistical phase calculated above. 

When the loops are knotted and linked, the exchange 
of two loops brings topological transformation of config- 
uration of the links. We first untie, unentangle the two 
loops to the simplest circle, and then exchange them. Fol- 
lowing an inverse braiding procedure, we restore the link 
configuration at their new position, so that one string 
evolves exactly into the same configuration of another. 
The Reidmeister moves are the fundamental propagators 
from one state to another in the Hilbert space of degen- 
erate ground state. In the three dimensional lattice, the 
exchange of two knots follows recursive Skein relation 
step by step. The statistical phase factor is no longer a 
simple complex number but a complex matrix inherited 
from the topological constrain. 

To operator on quantum states, we need to do braiding 
operations on the crossing point of Wilson loops one by 
one, and deduce the statistical matrix by using Skein re- 
lation. There is a different way to understand the mathe- 
matical braiding from physicist's point of view. As shown 
by Witten, the correlation function of Wilson loop op- 
erators in Chern-Simons theory is a topological invari- 
ant of Jones polynomial. We take two neighboring loops 
Lq = {Con Cfs) as one link configuration, and do a surgery 
by cutting one crossing point. Then we braid the inter- 
secting point on the cross section and glue them together 
again, it finally leads to new links, Li and L2- The parti- 
tion function of the new links and the original link follows 
the recursion relation AZ{Lo) + BZ{Li) + CZ{L2) = 0. 
We repeat this braiding operation on all the crossing 
points until all of the links are untied, then we exchange 
them, and braiding them up following an exact reverse 
operation. Finally we derive the algebra relation between 
statistical matrices. This is the whole strategy to obtain 
the statistical matrix algebra of links. 

The operations above sound awfully sophisticate, but 
it is applicable for both the continuous loop and dis- 
crete loops. As for the discrete loops, we have another 
different way to their algebraic statistical relation. In- 
spired by Levin and Wen's algebraic statistics relation for 
point particles hopping on two dimensional lattice [s^. 
we derived a general algebraic statistic relation of dis- 
crete loops using the projection operator Tf3a = | V'/j) (V'a | , 
which indicates the hopping from loop Ca to Cf}. If we ig- 
nore the orientation of the loop and take the flux quanta 
surrounded by the loop as one giant particles, the alge- 
braic statistic relation reads 

TapT^aTaX — TaxTjaTaf3 ■ (46) 

The statistic matrix strongly relies on the orientation 



of loops and the total number of particles confined along 
the loops. The total number of particles is not conserved 
when the loops transformed into another loop configura- 
tion. The hopping operator Tap must be represented by 
a matrix encoding the information of both the number 
of particles and orientation. Suppose loop Ca consists of 
n particles, and Cp consists of m particles, the hopping 
operator between Ca and Cp is 



/ hi 



ct Ci 



T, 



aP 



/i2i4q h 



hijc[c 



hikclck \ 

h2kC2Ck 



(47) 



\ h^iC^Ci hjijC^^Cj • • • hjikC^Ck / 



where c\, ■ ■ ■ , cj, are the creation operators of particles 
along loop Ca and Cj (j = 1, • • • , m) are the annihilation 
operators of particles along loop Cp . In three dimensional 
spin models, these fermions are introduced by three di- 
mensional Jordan- Wigner transformation Eq. (j29p which 
maps spins to fermion operators. The order of indices 
{i,j,...,k) encodes the orientation of Cp. For a special 
case that Ca and Cp have the same particle numbers, 
the hopping matrix is diagonal for loops with the same 
orientation and is skew-diagonal for loops with opposite 
orientation. The explicit hopping matrix on a torus of 
world sheet is 



/ hut 



lilll 





V 











\ 





(48) 



while the hopping matrix on a Klein bottle take the form 
of 



T, 



Klein 



( 











hlntln \ 



; 



(49) 



where iij = c[cj is the local hopping operator. Since 
Wilson loops project out many degenerate ground states, 
even if the loops are not entangled or knotted, the hop- 
ping operator does not always commute with each other. 



[fap,f^x]^0. 



(50) 



So there is no longer such a simple statistic algebraic 
relation as that for point particles in Ref. [32]. We have 
to develop new algebra on quantum braiding of knotted 
and linked loops, this exploration would be presented in 
some works elsewhere. 



VI. CONCLUSION 

In summary, a class of three dimensional spin liquids 
have been constructed and chiral spin liquids can be ob- 
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tained by following brief doping rules on conventional 
three dimensional lattices. For example, we introduced 
Kitaev-type 's coupling on the Cu-sublattice of crystal 
green dioptase Cu^SieOT, ■ 6H2O, and derived an exactly 
solvable spin model. Its exact ground is spin liquid. If 
we dope spin on its bonds and assume that there ex- 
ist only short range interactions between nearest neigh- 
boring doping spins, an exact chiral spin liquid emerges. 
The time reversal symmetry is spontaneously broken at 
ground state. Linked, knotted loop excitations are inter- 
esting elementary excitations obeying non-abelian statis- 
tics. The statistical phase factor strongly relies on the 
chirality of loops as well as the total number of parti- 
cles confined along the loop. Two loops with the same 
chirality behave like bosons. Two loops with opposite 
chirality obey fermionic statistics. Nontrivial multi-loop 



correlation measures topological quantum entanglement 
which is valid under continuous configuration transfor- 
mation. Any change of topology of links make the Link- 
ing number jump from one integer to another so that the 
topological correlation function correspondingly change. 
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